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Abstract 



This thesis contains two results for the low temperature behavior of quantum spin 
systems. First, we present a lower bound for the spin-1 XXZ chain in finite volumes 
in terms of the gap of the two-site Hamiltonian. The estimate is derived by a method 
developed by Nachtergaele in [22] called the Martingale Method. Our bound relies on 
an assumption which we have, as yet, been unable to verify analytically in all cases. 
We present numerical evidence that strongly indicates our assumption is valid. 

The second result is a proof that the spin- 1/2 d- dimensional XY model in the 
presence of an external magnetic field does not undergo a phase transition at low 
temperature, provided that the strength of the field is great enough. Using a contour 
expansion inspired by Kennedy in ^2j ; we show that the weights of contours satisfy 
a condition of Kotecky and Preiss ^H] which allows us to express the free energy of 
the system as a cluster expansion. As part of the setup we give a simple proof that 
the all-spin-up state is the unique ground state when the external magnetic field has 
strength at least 2d. 
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Chapter 1 



Introduction 



This thesis contains two results in the statistical mechanics of quantum spin systems. 
The first is an explicit lower bound bound for the spectral gap of the one dimensional, 
spin-1, XXZ ferromagnet. The estimate is obtained using a method developed by 
Nachtergaele in [22] • The second result is a proof on the absence of a phase transition 
for the Heisenberg XY model in an external magnetic field. In general the XXZ and 
XY models have very different behavior, ironically however, the techniques used in 
Chapter El were developed by Kennedy in ^2] to treat the low temperature XXZ 
model. 

The XXZ Hamiltonian on a domain A C Z d is given by 



where the operators are the spin operators in the direction e a at the site x G 
A, A > 1, and | • | is the standard £ 1 -norm on Z d . Ferromagnets are of general 
physical interest, and are not uncommon in nature. However there are other reasons 
to study ferromagnets as opposed to antiferromagnets or ferrimagnets. The first, and 
perhaps least serious, is that it tends to be very difficult to extract information from 
antiferromagnetic models. There are good results on antiferromagnets though. For 




(1.0.1) 



\x-y\=l 
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example in [7] Dyson, Lieb, and Simon used reflection positivity to prove the existence 
of long range order in dimension two or greater. However, the XXZ model is not, in 
general, reflection positive |2*5] . 

For the ferromagnet much more is known. In Alcaraz, Salinas, and Wreszinski 
gave a formula for the ground states of the XXZ chain for any spin. Koma and 
Nachtergaele showed that in the thermodynamic limit, the spin-1/2 XXZ model has 
a spectral gap above the ground state, and were able to calculate an exact expression 
for the gap [T^]. Later, Matsui used the existence of the gap to prove that for any 
anisotropy A > 1, the list of ground states given in £Q is indeed a complete list. 
Matsui's result was then improved upon by Koma and Nachtergaele in ^B] where 
they used the intrinsic quantum group symmetry of the XXZ chain to extend Matsui's 
result to the case A = 1. 

The existence of a quantum group symmetry of the finite volume spin-1/2 chain, 
initially noted by Pasquier and Saleur in j2H|, is still another reason for studying 
the XXZ model. Physcists generally associate any symmetry with physical meaning, 
while mathematicians tend to view it as a special but interesting situation that is 
deserving of further investigation. In ^3] Koma and Nachtergaele not only calculated 
the spectral gap of the spin-1/2 infinite volume chain, but proved a sort of ordering 
of energy levels (a la Lieb and Mattis in [T^]) for the finite chain. Specifically they 
showed that not only is the ground state of H* xz the L + 1 dimensional irreducible 
representation of SU q (2), but any first excited state must lie in an L — 1 dimensional 
irreducible representation. The quantum group symmetry does have some limitations 
however. It only exists in a restricted setting, i.e. for spin-1/2 chains in one dimension. 
In this paper we are interested in the spin-1 chain, so other tools must be used. 

The low-energy spectrum of the XXZ chain is also fairly well understood. This 
is relevant because the low-energy spectrum, and in particular the spectral gap, have 
direct implications for the low-temperature behavior, e.g. the decay of correlation 
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functions, of the physical systems modeled by the Hamiltonian. As stated before, 
Koma and Nachtergaele calculated the spectral gap of the spin- 1/2 chain in |15j . 
For spins greater than 1/2 Koma, Nachtergaele, and Starr showed that there is a 
non- vanishing gap in one dimension in [Hj. These results were extended by Caputo 
and Martinelli. In [Sj they showed that the gap grows linearly with the spin. This 
confirmed a conjecture of Starr in J3U]. The techniques of Caputo and Martinelli are 
limited however, in the sense that they do not give an explicit value for a lower bound 
of the spectral gap. In Chapter |2] we seek to remedy this in a small way. We present 
a method by Nachtergaele for estimating the gap of finite chains in terms of the gap 
of the nearest-neighbor interaction j^S^Sy + j^S^Sy + S^Sy. The method relies on 
an assumption (J2.4.8|) that we believe to be satisfied for the spin-1 (and larger J) 
chain. We have much numerical evidence to suggest that the assumption is satisfied, 
but not an analytic proof of the fact. For this reason the lower bound for the gap of 
the spin-1 chain is only a conjecture, but we are confident that it is accurate and will 
soon be rigorously verified. 

As stated before the second result in this thesis concerns the isotropic XY model 
in the prescence of an extermal magnetic field for dimensions 2 or greater. The 
Hamiltonian is given by 

{tc,a}eA x-gA 
\x-y\=l 

where the a a , a = 1, 2, 3 are the Pauli spin matrices and A is a finite subset of Z d 
with d > 2. In a paper dated 1961, [2D] , Lieb, Shultz, and Mattis studied the general 
anisotropic one dimensional, spin- 1/2, XY model without an external magnetic field 
(h = 0.) They rigorously established the equivalence of the XY model with a free 
fermion model via the Jordan- Wigner transformation, and proved the existence of 
long-range order for the anisotropic model using a method they developed for calcu- 
lating two point correlation functions. In the isotropic case the order vanishes. They 
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were also able to establish the uniquess of a gapless ground state in the thermody- 
namic limit for the isotropic case. Later, Dyson, Lieb, and Simon, using reflection 
positivity, proved that the spin-1/2 XY model has a phase transition in three and 
higher dimensions Finally, in 1988, Kennedy, Lieb, and Shastry demonstrated 
the existence of long-range order for the ferromagnetic XY model in dimensions two 
and greater, regardless of spin. This then proved the existence of Bose-Einstein con- 
densation for the hardcore gas on the same lattice at half filling ^3] • 

In Chapter |H] we show that for magnetic fields with h > 2d, there exists a range 
of non-zero temperatures where the XY model does not have a phase transition in 
two and higher dimensions. This case is not treated in the above papers since the 
prescence of the magnetic field destroys the reflection positivity of the model. In 
[T%] Kennedy and Tasaki study a much larger class of Hamiltonians, using similar 
techniques and with similar results. At this point, though, it remains unclear if there 
results cover the case with which we are interested. We employ a contour expansion, 
inspired by Kennedy in jT2], to show that the free energy of the system is analytic 
for some values of inverse temperature f3 and magnetic field h > 2d. Specifically we 
are able to show that the weights of our contours satisfy a condition of Kotecky and 
Preiss ^H] which allows us to write the free energy using a cluster expansion. 

I have been supported in part by the National Science Foundation under Grant 
No. DMS0070774, and VIGRE Grant No. DMS-0135345. 
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Chapter 2 

XXZ Model 

For many spin systems interesting behaviors (e.g. phase transitions) only occur in di- 
mensions greater than one. However, in the case of the XXZ model, phase transitions 
do exist in dimension one. This remarkable fact is closely connected to the symmetry 
of the spin-1/2 chain under the action of the quantum group SU q (2). 

In PH], Koma and Nachtergaele calculated the exact spectral gap for the spin-1/2 
XXZ chain. The goal of this chapter is a lower bound for the spectral gap of the 1- 
dimensional spin-1 model. This is obtained using a method developed by Nachtergaele 
in (221; which was later improved by Spitzer and Starr in [22]; for estimating the 
spectral gap of frustration free systems. 

We will start by introducing the Hamiltonian for the XXZ chain. Since the method 
we use to find a lower bound requires one to know the ground states of the finite 
volume Hamiltonian, we will introduce these, and prove that they are indeed ground 
states for all values of spin. For completeness, we also present a brief account of the 
infinite volume ground states, although this is a much more sophisticated subject. 
Next we will give a proof, due to Koma and Nachtergaele, that the spectral gap of 
the infinite volume spin-1/2 chain is given by 1 — A -1 , where A is the anisotropy 
parameter in the XXZ Hamiltonian. In the last section we introduce the method 
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used to obtain the lower bound and present our results. 
2.1 The Model 

Initially we consider only finite intervals A C Z. For each x G A we let 7i x = 
C( 2J+1 > be the configuration space for the individual spin at site x. The Hilbert 
space for the entire chain is given by H\ = (S^eA^r- The spin arises by way 
of a (2J+l)-dimensional irreducible representation of SU(2) on each 7i x . This is 
generated by three operators S%, S%, which are determined relative to a basis 
{\-J) x , |-J +!)„..., \J) x }by 

Sl\m) x = m\m) x (2.1.1) 

and 

S x \m) x = y/j(J + 1) - m(m ±T)|m± 1)^ (2.1.2) 
where S^r = S^. ± zS 1 ^. The 1-dimensional spin- J XXZ Hamiltonian is then given by 

jjxxz _ _ ('S'x'^+i + Sx^x+i) + •S'x'^+i (2.1.3) 

x=l 

where A > 1 is the anisotropy parameter. The overall negative sign makes this a 
ferromagnetic model, while the inclusion of the anisotropy forces low energy states 
to have large regions in which spins are aligned along the 3-axis. Notice that in the 
limit A — > 1 we obtain the isotropic Heisenberg, or XXX, model, while as A — ► oo 
the model tends to an Ising like Hamiltonian. 

As with many things in mathematics, the truly interesting behavior comes when 
we add boundary terms. To this end, we consider the so called "kink" boundary 
conditions. Let A (A) = \/T— A~ 2 and define H£~ by 

H+- = Hf xz + JA(A)(S 3 L - S'l) + J 2 (L- 1) (2.1.4) 
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The last term, J 2 (L — 1), is added so that the ground state energy is zero. If we let 
h x , x+x = -i (.S'.s;;.., + S^+i) + S^-i + JA(A)(Sl +1 + Si) + J 2 (2.1.5) 
then H~£~ can be expressed as the sum over translates of a single two-body interaction: 

L-l 
x=l 

The choice of A(A) is such that the spin-1/2 Hamiltonian will commute with 
representations of the quantum group SU q {2) where 2A = q + q~ l and q G (0, 1). If 
we define the operator K to be 

K = q- 2S3 (2.1.7) 
then K, S + , S~ satisfy the commutation relations 

[S + ,S-] = K ~ K 1 and KS ± = q^ 2 S ± K (2.1.8) 
q" 1 - q 

K and S ± are the generators of an irreducible representation of SU q (2). Amazingly, 
there is a way to define tensor products of representations of SU q {2). For a chain of 
L sites we write 



L 



S l,L] = lj ^-1 ® S l ® W ® ""® ^ ( 2 - L9 ) 

L 

5+ L] = ^ <g) • ■ ■ ® <g) ® l^+i <g> ■ • • ® Hi (2.1.10) 

2 = 1 

L 

^.i] = XI 11 ® " ' ® lx ~ x ® ^ ® ^"^+1 ® ■ ® (2.1.11) 

For spin-1/2 these operators commute with the XXZ Hamiltonian and were at 
the heart of the initial discovery of the ground states of the XXZ model. We will 
not use the quantum group symmetry much here and only mention the existence of 
the symmetry for completeness. Someone who is interested in learning more about 
quantum groups and there representations should consult the reference [TT] . 
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2.2 Finite Volume Ground States 

Perhaps the first, or most basic question, one can ask when given a Hamiltonian 
like ()2.1.4|) is: does it have any non-trivial ground states and if so, what are they? 
This question was first answered by Alcaraz, Salinas and Wreszinski in . As earlier 
indicated, they were able to realize the ground state of (j2.1.4|) in the spin-1/2 case 
as the highest weight irreducible representation of the quantum group SU q {2). They 
were able to generalize the result for spins J > 1/2. Here we present their result, but 
the proof provided is due to Starr in [3U] . 

Theorem 2.2.1 (Alcaraz, Salinas, and Wreszinski) For each eigenvalue m of^2 x=1 S x , 
m G { — JL, . . . , +JL} there is a unique (up to normalisation) state in the sector 
of magnetization m, such that H^ty™ = 0. ty™ is given by 

L ( 1 J \ l l 2 

*o m = E U ( i' x{J ~ mx \j +m ) \{™*}) (2- 2 - 1 ) 

{■m x }e{-J,...,J} L x=1 
2 a- rn x =m 

where \{m x }) is just the simple tensor 

\{ m x}} = \mi) ® \m 2 ) <8> ■ • • ® \m L ) 

with respect to the natural basis for the spin operators given in \2.1.1)) . 

Proof: The goal is to calculate all states (p such that H^~tp = 0, that is we want to 
calculate ker But by (|2.1.6|h is the sum of non-negative operators, thus 

if Lp G ker then ip E ker (h XtX +i) for 1 < x < L — 1. This gives 

L-l 

ker(if+-) = P| ker (^+0 (2.2.2) 

x=l 

Models with this property are often referred to as frustration free, since any state 
that minimizes the energy of the global interaction also minizes the energy of every 
local interaction and hence are not "frustrated" . 
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We begin by calculating f^^ker (h X)X+ i) for J = 1/2. For other values of J we 
embed the spin- J chain into a longer spin-1/2 chain and use representation theory to 
obtain the desired result. For J = 1/2, "Hl is the 2 L dimensional space ^ =1 C X . The 
simplest basis for TCl is given by 

L 

\{m x }) = \m x ) x 

x=l 

where {m x } £ {±1/2} L . Often we will use the convention that |l/2) = |f) and 
| — 1/2) = ||). It is a short calculation to show that 

h x , x+ i = -— (SlSl +1 + S%Sl +1 ) + SlS' x+1 + A(A)(Sl +1 - S^) + J 2 

is the orthogonal projection onto the the vector 

\0 = 7F^== {q\U)x,x+i - IIT)x,, + i) (2.2.3) 
V 1 + T 

where q is the unique solution to the equation 2 A = q + g^ 1 on the interval (0, 1). 
Thus, we can write 

^+i = IO(£l (2-2.4) 
To motivate our argument for general g consider the case when q — 1. Then 
|£)<b,x+i i s the spin singlet, or the anti-symmetric tensor, in TC x <S>H x +i. If r(x, x + 1) £ 
<Si is just the transposition of x and x + 1 then again for q = 1 we can write h x>x+ i 
as 

h x>x +l — — (l — -Pr(a),2i+1)) (2.2.5) 

where for 7r £ <S^, represents the standard action of <Sl on Ti L given by 

L L 

P.(g)|0x)x = (g)|^- 1 (,))x (2.2.6) 

x=l x=l 

Hence, for the isotropic model (A = 1), if ip is a ground state of h x>x+ i, then P T ( XtX +i)il} 
is also a ground state for all 1 < x < (L — 1). Since the nearest-neighbor transposi- 
tions generate the entire group Sl, we see that kernel of the isotropic Hamiltonian is 
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invariant under permutations, i.e. it is made up of symmetric tensors. This is quite 
evident when it is noted that the XXX Hamiltonian has full SU(2) symmetry, so 
that its ground state coincides with the heighest weight irreducible representation of 
SU(2) in ® L X=1 C 2 . 

For q 7^ 1 we again realize the ground states as symmetric tensors, however not 
with respect to the standard action P n . Let 
lo : {±1/2} L — > Z be given by 

L 
x=l 

and consider the basis for Hl given by 

\{m x })M=q»« m *»\{m x }) 

Now, let P( q ) : Sl — » GL(H L ) be the standard representation of the symmetric group, 
but with respect to the basis \m x )( q h Thus we have 

P^{m x })^ = \{m n - Hx) })^ 

but with respect to the standard basis we get 

pM\{m x }) = PM (q-»« m *»\{m x })W) 
= q-^ {m * ]) \{m^ {x) })^ 

This formula looks a little complicated, but it is exactly what we want. If we restrict 
our attention to the space T~C X (g> T~C x+ i then the symmetric group only consists of two 
elements, the identity e and the transpostition t(x,x + 1) = r. A short calculation 
shows that 

P (< ? )| TT) = | TT) p( 9 )| n) = g | iT) 
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Then we can immediately write down the symmetric tensors as 

ITT), |U>,iI±£# 

while the anti-symmetric tensor is given by 

IT1> — g|iT> 

Thus, by (I2.2.4|) . h x>x+ \ is the orthogonal projection onto the anti-symmetric tensor 
of the representation P^ q \ and hence we write 

h X)X +i = - — -Pt(2,z+i)) 

where the anisotropy parameter A in h x>x+ \ and q G (0, 1) are related by 2A = 
q + q' 1 . So, just as in the isotropic case, we are able to conclude that any frustration 
free state of H~£ ~ is invariant under nearest-neighbor transpositions. But nearest- 
neighbor transpostions generate the symmetric group, so the ground states of 
are symmetric with respect to the representation P^ q \ Since is identical to the 
standard representation with respect to the basis \{m x })^ we can immediately write 
down the symmetric states with respect to this basis: 



{m x }e{±l/2}L 

E rn x =m 

where m G {— L/2, (— L/2 + 1), . . . , L/2}. Expressing this in terms of the standard 
basis gives 



k = E (f {{mx}) \{ m *}) 

E m x =m 

= E <f- x{1/2 - mx) \{m,}) 

{m x }d{±l/2}L 

E rn x =m 

L 

E X[(f x{1,2 ~ mx) \{m x }) (2.2.7) 

{m x }e{±l/2} L x=l 
E rn x =m 
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which is the correct formula for J = 1/2. 

For J > 1 /2 we let H { X J) = C 2J+1 as above. H { X J) can be realized as the highest- 
weight space in the representation CS^C^.. In other words Tix^ consists of the 



symmetric tensors in 
want is 



tf=iHx/ ■ But we have L copies of Ti^ thus the space we 



2.1 



(J) 



<g>(®«2? 

i=i \j=i 

Again, suppose that P : S 2 j — > GL(H^J 2 ^) is the standard representation of S 2 j 
on 'H^] 2 ^ with respect to the basis \ {m Xj }). Let S^j be the direct product of L copies 
of the symmetric group on 2J elements. Define a representation II of S£f onto 



(2.2.8) 



by 



^ u ...,7r h )Q$\{m Xj ,i})i 

i=l 



i=i 



so that the i th copy of S 2 j acts in the standard way on the i th copy of in 



(J) 



Let Sym n [IHI^] C H^ J ^ denote the subspace of M.^ which is invariant under 



the action of II. This space coincides with Ti.^ because 



Sym n [H 



(Jh 



Sym n 



2.1 



i=l \j=l 
2J 



(g) Sym, 



8=1 



.i=i 



(2.2.9) 
(2.2.10) 



since II only acts on individual factors in via P. But 

" 2J 

Sym P 



Xj,l 



U=i 



ft 



(•/) 



and hence 



Sym n [4 J) ]=(g)H 



(J) 



8=1 



Next, we define a spin-1/2 Hamiltonian on H^ J ^ which reduces to the Hamiltonian 



. We consider L distinct and disjoint copies of the set {1,2,..., 2 J} which we 
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denote by [1, 2 J]j, f = 1, . . . , L. Then define H^uf^2j L ) ^ 



L-l 

,(1/2) 



H (2jtl.lj L ) = Y1 E K,x i+ i (2.2.11) 



t— 1 (x it x i + 1 ) 

Xie[l,2J]i, x i+1 £[l,2J] i+1 

Hpjf 1 ^^ is clearly invariant under the action of II. Therefore, setting Psym n as the 
orthogonal projection onto Sym n [H^], gives 

TT+-AJ) _ P TT+-, (1/2) p 

L — r Sym n - n (2J 1 ,...,2J L ) r SD"'n 

. However, the construction of II was such that H^ 2J ^ 2 2J ^ commutes with the entire 



action, thus 



TT + -,(J) _ rH— ,(1/2) p 

n L — S1 (2J 1 ,...,2J L ) r Sym u 



This means that kerfif^™ ) = ker(if 2 l J 1 2 < ' 1 ^ 2 ^)riran(P sym u ) However, by the argument 
for the spin- 1/2 case, the ground state of 2J L ) * s invariant under the action 

P^ q \ Therefore 

M#J^!!5j l )) c ran(P 5j/mn ) 

which implies 

ker(tfr ,(J) ) = ker(iy+;; ) (1 / 2 \ ) ) n ran(P Syfnn ) = ^{H^% l) ) 

that is, the ground state of is just the ground state of Hnjf 1 2j L y But the 

formula for the ground states of Ht 2Ji 2 Jl \ is given by ()2.2.7|) . Writing these in terms 
of symmetric tensors; doing so yields 

2JL 

k = e n^ (i/2 ~ m * j) iw> 

{m :E;) }e{±l/2} 2JL J=l 



Ej m Xj =m 



2JL 



{m^}e{±l/2}2Ji J=l 
Ej m Xj =m 
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Next, we rewrite the single sum as two sums; first over the L chains of length 2 J, and 
then over the bases of the individual chains. This gives 

/ \ 



n 



{m x }e{-J,...J} L z=l 
Es m x =m 

L 

E 

{m x }e{-J,...J} L x=1 
Yj X m x =m 

e n< 

{m x }e{-J,...J} L x=l 

Y^x m x =m 



-x(J—m x ) 



-x(J—m x ) 



-x(J—m x ) 



E 



2.1 



{^€{±1/2} j = l 
\T-j b j,x=m x 



2J 

J + m. 



2J 

J + m x 



1/2 



|m x ) 3 



1/2 



which is the same as (|2.2.1j) . 



2.3 Infinite Volume Ground States 

In infinite volume the problem of determining the ground states for the Hamiltonian 
is much more subtle; even defining the Hamiltonian turns out not to be as simple as 
one might think. For AcZ, |A| < oo the algebra of observables is A\ = B(H,a). In 
infinite volume the algebra of quasi-local observables Aoo is defined by 



A 



U ^ 

ACZ 

\|A|<oo 

where the completion is taken in the operator norm topology. Note that if A' C A 
then there is a natural way to embed A\' as a subalgebra of A a, namely 



/ 



A A > 3 A' ^ A' ® ]1 A \a' € A 



(2.3.12) 



In the infinite volume H% is determined by a derivation 5 on Aoo. F° r an y A G Aa, 

|A| < oo 

6(A) = lim \H A i,A] 
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A state, uj, on Aoo is a positive (normalized) linear functional. A ground state is one 
that satisfies the local stability property 



This condition becomes much more intuitive if we go to the GNS representation 
corresponding to the state u>: (it, T~C, Q). For the XXZ Hamiltonian with H — boundary 
condition the GNS Hilbert space, H.+- , can be understood in a simple way. Let Q + ~ 
be given by 



is a dense subset of 7i + _. For a local observable A, the stability condition (J2.3.12)) 
becomes 



uj(A*5(A)) = (Q + ^\7t(A*)5(7t(A))Q + ~) 

= (Q + '\tt(A*))[tt(H),tt(A)}Q + ' 

= (Q + ~\7t(A*)7t(H)tt(A)Q + -) - (Q + -\tt(A*)tt(A)tt(H)Q + -} 

= (Q + ~\7t(A*)7t(H)tt(A)Q + -) 

= (tt(A)Q + ~\ tt(H) \tt(A)VL + ~) > 



Thus we see that (12.3.12)1 designates a ground state as one for which any local change 
coincides with some increase in energy; which is exactly what a ground state should 
be. 

A closely related, but stronger, notion in the infinite volume is that of a zero 
energy state. For a model like the XXZ, we can define this in the following way. 
Let H\ = J2{ x i+i}ca hx,x+i where h x>x +i is the translate of some nearest-neighbor 



uj(A*S(Aj) > A e A\; |A| < oo 




(2.3.13) 



Then the set of vectors of the form 




(2.3.14) 
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interaction for all (x,x + 1). Then u is a zero energy state iff 

u(K, x+ i) = min u)'{h X)X+1 ) 

for all {x, x + 1} C Z where here the states are all normalized. This is also called 
frustration free as in the finite volume case. If A e .A[i,l], then 

u(A*6(A)) = ]im u ( A* ( H A A- AH A )) 



^2u(A*(h x>x+1 A - Ah XjX+1 )) 

x=0 
L 

J2u(A*(h x>x+1 -\ )A) 

L 

^2^A(h x , x+ i - A ) > 



x=0 

where Ao is the smallest eigenvalue of h XtX +\ and uja{-) = w(A* ■ A). Therefore, if 
u is a zero energy state then it is also a ground state, however the reverse is not in 
general true. For the XXZ chain Gottstein and Werner found all zero energy states 
in [H]. Their results and techniques were quite general and will not be presented 
here, although they obtained very nice results connecting zero energy states and VBS 
states. Moreover, they predicted correctly that for the XXZ and XXX chains the zero 
energy states are the only ground states. Recalling the definition of in (|2.3.13|) . 
the zero energy states for the XXZ Hamiltonian are given by 

oo k n+k 

^ (n) =J2 q-fr+-+**)+<K+-+v<*») Y[ S~ Y[ S+n + ~ (2.3.15) 

k=0 {xi,...,x fe }CZ< J = 1 i = 1 

{!/li-i!/n+)t}C2>l 



It was shown by Matsui in [21] that the zero energy states defined above along with 

those obtained from — h, ++ and boundary conditions account for all ground 

states of the XXZ chain. Later, Koma and Nachtergaele in JH] gave an alternative 
proof of this fact with the added result that their techniques could be applied to the 
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XXX chain as well. This is a major difference between Koma and Nachtergaele's and 
Matsui's results. Their results not only apply to spin- 1/2, but to all J. 

2.4 The Spectral Gap 
2.4.1 Spin-1/2 

The spectral gap of the spin-1/2 XXZ chain was first calculated by Koma and Nachter- 
gaele in ^5] • Their proof is a truly beautiful piece of work that makes use of standard 
tools of statistical mechanics, such as the transfer matrix, while also employing newer 
phenomena such as the quantum group symmetry of the XXZ model. Together with 
S. Starr they later showed for all J there is a non-zero spectral gap in the thermo- 
dynamic limit, although they did not provide an estimate ^7j. For completeness, we 
provide the proof of the existence of the gap in the spin-1/2 case. 

For a, [3 £ {+, — }, let Vl al3 be the generalizations of Q + ~ given in ([2.3. 13|) . that is 



la) if x < 

fi a ^ = (S?)n a * 3 (x) where = { (2.4.1) 



if x > 1 

and similiarly the GNS space 7i Q/ 3 has a dense set of vectors of the form 



V> = AQ af3 where A £ (J A A , |A| < oo (2.4.2) 



The cyclic vector of the GNS representation, £Iq NS is then given by a generalization 
of equation ()2.3.15|) 



nf NS = Z~\q) E 9 E ^ (W -^ II S^tt* (2.4.3) 

fc=0 {x lr ..,x fc }a< j=l 

fel,-,!/!:}^! 

where Z~ x (g) is a normalization factor. The Hamiltonian H^f NS is then defined by 
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the equation 

H&gAn&a = }f^ H f' A ^gns (2-4.4) 
Note that the set of vectors in (J2.4.2|) are in the domain of Hqns- 

Theorem 2.4.1 (Koma and Nachtergaele) For all A > 1 and for any choice of 
GNS sector ({+-), (-+), (++), ( )) let T af3 be the set of all A G M such that 

(Q al3 \A* (HgnsY AQ a(S ) > \(Q al3 \A* (h^ n X AQ a(3 ) (2.4.5) 

for all A G Ua^4a- Then 

sup r a ^ = 1 — A 1 (2.4.6) 
or, in other words, if 7 Q/ g is the spectral gap of H^f NS then 

la p = 1 - A" 1 . (2.4.7) 

The proof of theorem 12.4.11 consists of two main parts: first an estimate on a lower 
bound for 7 and then an argument to show that this is also an upper bound. 

Proof: In order to bound 7 from below by 1 — A -1 we show for all local observables 
A 

{o. aP \A* (Hf Ns y An*?) > (1 - a-^^ia* fe^ s ) 2 AQ a/3 ) (2.4.8) 

that is, for every ip in the range of H^f NS 

(^^^^(l-A- 1 )^!^) (2.4.9) 
But, if A G ^4a with A finite, then 

(tt#\A* (Hf NS ) k Att#) = (SV#\A* (H^ im ) k A^) (2.4.10) 
The operator A* fi^^-i ^ j A is in ^4a+{-i,o,i}' an d thus we can write 

(U#\A> (Ht^_ 1M ) h AO*) = Tr L +{ -i,o,i}^ (^+{-1,0,1})^! ( 2 - 4 - n ) 
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where Pa+{-i,o,i} is a density matrix. Hence, we have shown the infinite volume gap 
7 is greater than or equal to 1 — A^ 1 if, for all finite volumes A, the spectral gap of 
H^f is greater than 1 — A -1 . This is the result of the following proposition. 

Proposition 2.4.2 Let denote the spectral gap of H£~ ■ Then for all A > 1 and 

L > 2 we have 

7i = 1 - A" 1 cos(tt/L) (2.4.12) 

To prove Proposition 12.4.21 we calculate the first non-zero eigenvalue in the subspace 
of TIl corresponding to eigenvectors of S'^ ot with eigenvalue L/2 — 1, i.e. the sector 
with one downspin, and show that this is also a lower bound for the smallest non-zero 
eigenvalue in all other sectors. 

Let D x = S~\1) L where ||) L = ®^ =1 ||) x . Then if ^ is a vector in the sector with 
one downspin 

L 

i> = ^a x D x (2.4.13) 

x=l 

Now suppose that ip is an eigenvector with eigenvalue 8, then 

L 

H+'xj; = a x H+~D x = 8^ (2.4.14) 

x=l 

On the other hand computing H£~ gives 

„ t - Dl = LiM A -1- D2 

L 2 2A 

H+-D x = D X -^(D X ^ + D X+1 ) ioi2<x<L 

H + L -D L = 1 -^^-n L -±D L _ 1 (2.4.15) 

Putting (|2.4.13|) - (j2.4.15|) together we get the equations 

a x+ i = 2A(1 - E)a x - a x - X for 2 < x < L - 1, (2.4.16) 
a 2 = 2A(l/2 + A(A)/2-£)a 1 (2.4.17) 
a L _i = 2A(l/2-A(A)/2-£)a L (2.4.18) 
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Notice that if ax = then the above implies tjj = 0, so we can assume from now on 
that a\ 7^ 0. Equation (J2.4.16)) can be written as 



(2.4.19) 



2A(l-£) -1 

where T — | | is called the transfer matrix. Iterating this gives 

1 



/ a x+ i i 




' ax \ 




M 




\ a x ) 









| = T L ~ 2 1 








I a, / 



(2.4.20) 



So taking into account (j2.4.17jl . <|'2.4.1H|) . and (j2.4.2()jl we arrive at 

J 1 U^H 1 *^ ,,4,!) 

\ A[l - A(A) - E] J \ 1 / 

Equation (j2.4.21|) can be solved using the eigenvectors and eigenvalues of T. The 
characteristic equation of T is 

A 2 - 2A(1 - 8)\ + 1 = (2.4.22) 

which has roots 



A± = A(l-£)± VA 2 (l-^) 2 -l (2.4.23) 

Initially, let us consider the case when T does not have degenerate eigenvalues, that 
is A(l — 8) 7^ ±1. The eigenvectors of T are 

e± = ( X * j (2.4.24) 

We want to solve equation ()2.4.21|) so we must decompose each vector into a sum of 
e + and e_. Projecting the vector on the left of (j2.4.21|) onto e± gives 

1 



A[l - A (A) - 8] 



/i + A_e + + /i_A+e_ (2.4.25) 
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where 



^=2 



1 ± A- 



A(A) + £ 



'A 2 (l-£) 2 -l 
and the one on the right hand side of (j2.4.21|) gives 



(2.4.26) 



A[1 + A(A) -£} 
1 



r] + e + + ?7_e_ 



where 



1 ± A- 



A(A) — £ 



^A 2 (\-£) 2 - 1 

Now we can substitute 1)2.4.25)1 and 1)2.4.27)1 back into ()2.4.21)1 to get 



a^(/i + A_e + + /i_A + e_ 



a 1 T L_2 (?7 + e + + ?7_e_) 
a 1 (r/+T L - 2 e+ + r/„T i - 2 e_ 



But e + and e_ are linearly independent thus we get the equations 

,L-2 



(2.4.27) 
(2.4.28) 

(2.4.29) 
(2.4.30) 
(2.4.31) 



ai 



A 



— A 4 - 
ai 



A_ 



A. 



A- 1 
A i_1 



(2.4.32) 
(2.4.33) 



where we have used the fact that A + A_ = 1. Simplifying (J2.4.32)) and ()2.4.33|) further 
gives 

(2.4.34) 



x 2L-2 = V± X ± 



where we have assumed that r/_ 7^ 0. This assumption is not restrictive since 



?7_ = 



1- A- 



A(A) — £ 



'A 2 (l -£) 2 - 1 







1 

2 

A[A(A)-£] = ^A 2 (l-£) 2 -l 



£VA 2 -1=A£ 
£ = 



(2.4.35) 

(2.4.36) 
(2.4.37) 
(2.4.38) 
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thus r]_ = corresponds to the ground state in the L/2 — 1 sector. Moreover, r]_ ^ 
guarantees that all /i± and r]± are non-zero, so we see that equation (|2.4.34|) is indeed 
valid in the setting we want. Using (J2.4.26j) . (|2.4.28|h and (J2.4.23j) . we can write 



, + _A + v ^T x A + -(A-7A^T) 



A*- A+ A+ - (A + VA 2 -1 

and 

A - VA 2 - 1 A+ - (A + VA 2 - 1 



(2.4.40) 



V+ A+ A + — (A - v/A 2 ^!) 

Together with ()2.4.34j) this gives 

Xf = 1 (2.4.41) 

which implies 

A + = e M l L £eZ (2.4.42) 

Returning now to ()2.4.23j) and solving for Sl(£) gives 

£ L (£) = 1 - A+ 2 "^ A " = 1 - A" 1 cos(tt£/L) £ = 1, 2, . . . , L - 1 (2.4.43) 

where the condition on £ comes from the assumption that ()2.4.22|) has non-degenerate 
roots. Furthermore, we see that equation (|2.4.43|) gives L — 1 distinct eigenvalues, 
which gives us a total of L distinct eigenvalues of in the (L/2 — 1) sector when 
we add the ground state energy £ = 0. More to the point, there are not any solutions 
with A(l -£) = ±1. 

The above proves Proposition l2.4.2l that is in the (L/2—1) sector the smallest non- 
zero eigenvalue is given by 1 — A -1 cos(7r/L). To finish the proof of Theorem 12 . 4 . 1 1 we 
must show that this is same for all other sectors. In order to accomplish this we utilize 
the quantum group symmetry of the model. Note that any irreducible representation 
of SU q (2) of dimension L — 1 must intersect the sector with one overturned spin. Thus 
if the first excited state, (p™, lies in an L — 1 irreducible representation for any choice 
of magnetization m, 1 < m < L, we can raise it to a state, 0™, in the L/2 — 1 sector by 
way of the quantum group raising operator. The state cj>™ has energy greater than or 



CHAPTER 2. XXZ Model 



23 



equal to 1 — A -1 cos(7r/L) by the argument above. However, since the quantum group 
commutes with the Hamiltonian H~£~, <f>™ has the same energy as the original ground 
state, <p™, and Theorem 12.4.11 follows . Thus we must show that the lowest excitations 
lie in the L — 1 dimensional irreducible representations of SU q (2). To this end we 
present the following lemma without proof. The proof makes use of straighforward 
applications of the quantum group symmetry as well techniques similar to those in 
Section 12X21 

Consider an arbitrary spin chain of length n. Let 7i S 3 t > n j_ s be the span of 
eigenvectors of Sf of with eigenvalue greater than nJ — s (less than nJ.) Let H n = 
X^=i hx,x+i be an SU q (2) invariant Hamiltonian on a spin J chain with n sites. Let 
7 n denote the spectral gap of H n and let en be given by 

£ W = min «^» (2.4.44) 



Then the following holds 



Lemma 2.4.3 (Koma and Nachtergaele [15J) Let Hl be as above. Suppose that 

h x ,x+\ > and that ker(if^) is non-empty. Furthermore, suppose that ker(iJ„) co- 
incides with the irreducible representation of SU q (2) with maximal spin (=nJ,)for 
2 < n < L. If 

eW < e^l (2.4.45) 

for all 2 < n < L — 1, then 

1L = e? J) (2.4.46) 



The calculation with the transfer matrix, combined with Lemma \2. 4. 31 completes the 
proof of Proposition 12.4.21 Thus, for all L > 2, the spectral gap of is given by 

7l = 1 - A" 1 cos(7r/L) > 1 - A" 1 (2.4.47) 
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Thus 1 — A^ 1 is indeed a lower bound for 7. 

In order to complete the proof of Theorem 12.4.11 we must show that 7 is less than 
or equal to 1 — A -1 . This is a little subtle; upon initial inspection one can see that 
the gaps for Hq^s anc ^ ^gns are the same by s P m flip symmetry, and that the gaps 
for Hq^s an d Hqns are the same by left-right symmetry. However, the ground state 
of Hq^ s is non-degenerate but the ground state of H^ s is infinitely degenerate, so 
it is not so clear that the gap for and the gap for Hq^ s are the same. It is 

known that the translation invariant states are weak limits of the (anti-) kink states. 
Thus, 

. (Sl°P\A%H£ { _ im ) 3 AQ*) 

a™a a (n^\A*(H+- { _ 1M ) 2 An^) 

- A,^A A «i±oo (WP\T n (A*{H+- { _ im ) 2 A)WP) 

(Q ++ \A*(H+7r lon ) 3 AQ ++ ) 
AeA A (0++|A*(F+- { _ 1A1} )MO++) 1 ' • ] 



A 



so 7 + ~ < 7 ++ and hence we must only show that 1 — A -1 is greater that or equal 
to 7 ++ in order to finish the proof of Theorem 12.4.11 Towards this end we start with 
the variational priciple 



7 ++ = inf ^'fffi 5 ^ (2-4-49) 



i/>Gdomii' ( t^ s 



In addition we note that since h x ^ x+ \ > 0, for any sequence of sets Ai C A2 C 
■ ■ ■ , n G N, H.++ D kerffjj" 4 " D kerif^" 2 + D • • ■ . So, if iJj must be in a space orthogonal 
to ker Hqnsi ^ nee d only satisfy ip _L keri?^ + for some particularly well chosen A. 
Let A = [1, n] C Z. Define the spin wave operators, A% fc G —{1, 2, ... ,11— 1} 

by 



l k 

= —J2e ikx S- (2.4.50) 



x=l 
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where we multiply by 1 / ^fn so that 

(n ++ \X* k X l n ++ ) =5 k)l (2.4.51) 

To get the upper bound we take ip = (c\Xk 1 + c 2 Xfc 2 )f2 ++ . Notice that ip is in the 
subspace of 7i++ with one overturned spin; specifically, Y22=i ^x^P = ( n /2 — l)^- But 
for each eigenvalue m of YlZ=i ^x there is exactly one ground state tp™ of Ht ~, such 
that Xl™=i ^x^cT = m i )1 Q- Therefore for fixed distinct k\ and fc 2 we can pick c\, c 2 
such that ip J_ ip^" 2 ^ 2 . But this implies that 

n,fci,fc 2 Cl)C2 ||^| 

However, we will show that 



inf sup w ™ 2 syJ/ > 7 ++ (2.4.52) 



inf sup = 1 — A -1 (2.4.53) 



which together with (|2. 4.52(1 establishes 1 — A 1 as an upper bound for 7 ++ . 
For 1 < x,y < n define T x <y by 

T x , y = (Q + \S x H^ n+1 ^S y f] ++ ) 

= ^(2A^ - fa, y - 1 - 5 XiV+1 ) (2.4.54) 

Then a simple, messy computation shows that the sup Cl C2 in (|2.4.53(1 gives the norm 
of the matrix M(n, ki, fc 2 ) whose matrix elements 
M(n, k\, k 2 )i,j = M n (ki, kj) are given by 

M n (k,l) = ^it e ~ tkXT *,y eilV = M 1 " A^ 1 cos A;) + (2.4.55) 

x,y=l 

for k, I = 2irm/n. Thus we have 

inf sup ^fffi/^ = inf ||M(n, fci, fc 2 || = 1 - A -1 (2.4.56) 

n,fci,fc 2 ci,C2 r n,fci,fca 



which finishes the proof that 7 ++ < 1 — A x , which in turn gives us Theorem 12.4.11 
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2.4.2 Spin-1 

For spins greater than one-half, the quantum group symmetry, which was instrumental 
in the calculation of the spectral gap in the spin-1/2 model, is lost. Although this 
makes questions about the gap for spins greater than one-half harder, it does not make 
them untreatable. In Koma, Nachtergaele, and Starr showed that for all values of 
spin, there is a non-vanishing gap in the thermodynamic limit. They also conjectured 
that the gap grew linearly with the spin. This prediction was confirmed in [5] by 
Caputo and Martinelli as a corollary to results on interacting particle systems. Using 
the equivalence of the Markov generators for certain react ion- diffusion processes and 
Hamiltonians of quantum spin models, Caputo and Martinelli were able to show that 
for all values of anisotropy parameter A > 1 there is a constant 5 such that 



for all values of spin- J and chains of length L > 2. Here we use techniques developed 
by Nachtergaele to give a concrete lower bound for the gap in the case J = 1. We 
begin by introducing the techniques for a slightly more general class of Hamiltonians 
than the XXZ; specifically the method will work for any Hamiltonian with frustration 
free ground states, so we only assume this. The method presented here has been 
significantly improved by Spitzer and Starr in but there method is a bit more 
complicated, so for simplicity we present the original method here. We then give the 
estimate and then seek to justify it through a series of Lemmas. 

The Martingale Method 

We start with a spin chain on L sites and label the Hilbert space as TLl = (C d )^ L . 
Let hip > be a two-site Hamiltonian acting non-trivially on the first two sites, 
and let h x>x+ i denote the translate of hx s 2 acting non-trivially on the factors at sites 
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{x, x + 1}. We consider the Hamiltonian 

L-l 

Hl = h XtX+ i 

x=l 

and assume that the kernel of Hi is non-trivial. \il<a<b<L then we denote by 
G[ a> b] the orthogonal projection onto the kernel of Y^ x ~= a hx,x+i, and use the convention 
that G{ x y = 1. More generally, for A C [1,L] we set G A to be the orthogonal 
projection onto 

ker ^ h x,x+i 

x,{x,x+l}cA 

Then the 67 a satisfy 

G Al G A2 = G A2 G Al =G A2 HA 1 cA 2 (2.4.57) 
G Al G M = G A2 G Al if AxnA 2 = (2.4.58) 

The above follows directly from the assumption that h 12 > 0, and thus we have for 
any A C [1,L] 

ker H A = ("} ker h x>x+ i 

x,{x,x+l}cA 

Moreover if 7 is the smallest non-zero eigenvalue of hi 2 = H 2 then 

h x , x+ i >7(l-G [a!|X+ i]) (2.4.59) 
Next define the family of orthogonal projections {E n : 1 < n < L} by 

H — ^[1,2] if n = 1 

G[i,n-i] - C?[i,„] if 2 < n < L - 1 (2.4.60) 
G[i, n ] if n = L 

Using (|2.4.57|) and (|2.4.58|) it is clear that the E' n s are self-adjoint and satisfy 

E n E m = 5 n<m E n (2.4.61) 

L 

Y,E n = 1 (2.4.62) 



n=l 
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Theorem 2.4.4 (The Martingale Method) Assume that H is a Hamiltonian sat- 
isfying the properties above. Furthermore, assume that for each n, 1 < n < L — 1, 



\\G [n>n+1] E n \\ < 1/V2 (2.4.63) 
Then for all if) such that G[\^ n \^ = 

(VW> > 7(1 - V2e) 2 \\iP\\ (2.4.64) 

where e = max n ||G[ nin+ i].E' n || and 7 is the spectral gap of the nearest-neighbor Hamil- 
tonian hi 2- 



Proof: Take tp such that G [ljn] i/j = 0. By (12.4.611) and TZTW, 

L 

ihi 2 = Eii^ii 2 

n=l 

but E n = G[ 1>n ] so 

IHI 2 = EII^II 2 

n=l 

We then use (ip\h XjX+ iip) to estimate H-En^ll by writing: 

II^H 2 = ('^1(1 - G [nin+1] )E n il>) + (4j\G [n , n+l] E n iij) 

then, we insert the resolution of the identity given by (j2.4.62j) into the second term 
noting that G[ nin+ i]ip = to get 

L-l 

\\E n ^f = - G [n , R+11 )^) + (^I^^G^+i,^) (2.4.65) 

m=l 

The second term simplifies further. 

E m G 

[n,n+l] (G [l,m] (J ' [l,m+l] 

— C[ n:n+ i](G?[i jfn ] — G?[i im+ i]) 
= G\ nn+ i]E m 
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for m < (n - 1) or m > n + 1 by (j2.4.58|) . But if m ^ n, then E m E n = by ()2.4.61() 
so that 

||£„V|| 2 = {VKl-Gln.n+l])^) 

+ {t(f\(E n .i + £ n )G [n , n+1] £ n ^) (2.4.66) 
= (^\{l-G [n;n+1] )E n ^) 

+ ((En-! + S n )^|G [n , n+1]J B n ^) (2.4.67) 

We now estimate both terms in (|2. 4.67(1 using the identity 

Kvil^>l<^IW 2 + ^IW a 

for c > 0. Hence 

\\E n ^\\ 2 = ((11 - G [ntn+1] )i/>\E n t/>) + ((£„_! + E n )ip\G [n!n+ !]E n ijj) 

< -L^KH-C^^ + I^l^) (2.4.68) 

+ — (^I^G^+i]^) + f (V'K^-i + K)V) 

= ^(VKl-C^+nW + lll^H 2 (2.4.69) 
+ ^||G [n , n+1] ^f + |||(£ n _i + E n )^l| 2 

The first term in ([2.4.69(1 . ^-(^(l — G[n,n+i\)' l l } ) > can be estimated using ((2.4.59(1 and 
gives 

^-(if>\(l - C [n!n+1] )^> < ^—(^\h n , n+ !^) (2.4.70) 
2ci 2ci7 

The third term is then treated using (|2.4.63|) to get 

i|G [n , n+1] i^|| 2 < ^\\E r M\ 2 (2-4.71) 

ZC2 /C2 

Finally the last term can be split up as 

f||(£n-i + K)V|| 2 = f (IlK-i^ll 2 + ll^f) (2-4.72) 
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since E n _i and E n are mutually orthogonal projections. Thus by (|2.4.69() - (j2.4.72() we 

have 

2ci7 2 ^ c 2 

+ |(||K-i^l| 2 +II^H 2 ) 



Summing over n gives 



1 / / | TT I \ , ( C l , ^ 2 



< — <^|^)+ - + — + c 2 
2ci7 V 2 2c 2 

where we note that there is not an Eq term and that Ejj\\) = 0. Rearranging this 

expression gives the inequality 

(iP\Hi;) > 7Cl ^2 - c x - ^ - 2c 2 ^) II^H 2 (2.4.73) 

Maximizing C\{2 — c\ — ^ — 2c 2 j over ci > 0, c 2 > gives c\ — 1 — ev^2 and 
c 2 = e/v^- Subsituting these into (12.4. 73j) gives (I2.4.(i4l) . 



Main Result and Conjectures 

We present here an explicit lower bound of the spin-1 XXZ Hamiltonian based on the 
Martingale Method given above. The result relies on an, as yet, unverified assumption 
()2.4.63|) . Since the projections G\ and E n all commute with S^ ot = ^2^=1 ^ii we can 
first project into each sector of S^ ot , and then check to see if assumption ()2.4.63() is 
true. This is the strategy we employ and it is the substance of Lemma 12.4.71 and 
conjecture (|2.4.8j) . We are able to verify that for a spin-1 chain of L + 1 sites that 
()2.4.63|) holds in the sector of magnetization L. For the remainder of the sectors 
we are unable to verify (|2. 4.63(1 analytically, but we do have numerical evidence to 
support it. 
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We proceed as follows: In Theorem 12.4.51 we apply the Martingale Method to 
get a lower bound for the gap of the spin-1 XXZ Hamiltonian as a function of the 
two-site gap 7. In Lemma [2.4.61 we calculate the gap, 7, of the two site Hamiltonian. 
In Lemma 12.4.71 we calculate ||G[l,l+i]-Ex|| in the sector with one overturned spin 
and show that it is less than l/v2 for all L, followed by the conjecture, along with 
numerical support, that the maximum ||G7[l,l+i]-Ex|| is achieved in the sector with 
one overturned spin. 

Theorem 2.4.5 Consider the spin-1 XXZ chain on L + 1 sites, and assume that 
Conjecture 2.4-8 holds. Let ip 7^ be such that G^l+i}^ — 0- Then 



where 7 is the gap of the two-site Hamiltonian given by | — y § — ^r- 
Lemma 2.4.6 (7: gap of ^1,2) Let hx s 2 be given by 

h,2 = ~ {SlSl + S* S 2 2 ) - SfSi - A(A) (Si - Si) + 11 
and let ip 7^ 0. Then 



i>:d{;;^=Q uf 2 V4 a 5 



7 = , min ^ ^IX?> = --\l--— (2.4.75) 



Proof: This is really more of a calculation than a proof. The calculation proceeds 
in the following way: Let P m be the projection onto the subspace of Ti.2 with total S 3 
component equal to m. We calculate 7 m = gap (P m ^i,2-fm) for m = —2, —1, . . . , 2, 
which is just the first non-zero eigenvalue since hx s 2 has a non-trivial kernel. Then 7 
is just the minimum of the 7 m . For m = ±2 the subspaces are one dimensional and 
both vectors are ground states. For m — ±1 both P±xhx,2P±x can be expressed as the 
matrix 

' 1 + A(A) -i 
P ±1 h h2 P ±1 =\ J (2.4.76) 

-i l-^(A) 
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where A(A) = a/1 — 1/A 2 . The matrix has eigenvalues and 2, so we see that 



7±i = 2 



(2.4.77) 



For m = 0, P^hi 2P0 the matrix is 



( 2(1 - A(A)) 



i 1 



V 



(2.4.78) 



-i 2(1 + A(A)) y 

Here the matrix also has a zero eigenvalue; when A is factored out of the characteristic 
equation we get 



A 2 -5A + (— +4) = 



This has roots 



so we get 



x 5 + 9 2 



5/9 2 
70= 2-Vl"A^ 

Taking the minimum of (|2.4.77|) and ()2.4.79|) we see that 



_ 5 _ 9 _ 2_ 
7 ~ 2 " V 4 ~ A 2 



(2.4.79) 



(2.4.80) 



We now show that in the sector with magnetization L, \\G[l.l+i]El\\ < l/v2 for 
all q G (0,1). This is the first step in verifying ()2.4.63|) for the spin-1 chain. We 
mention that the calculation given here for the spin-1 case immediately generalizes to 
higher J, but we do not present that here. Let P m be the orthogonal projection onto 
the sector with total S 3 component equal to m. Then define e m = \\P m G[L,L+i]ELPm\\ ■ 



Lemma 2.4.7 el 

(0, 1), e L < 1/V2. 



(l-g 2L ) q 2 
(l-q 2L + 2 ) ' 1+q 2 



,1/2 


r q 2 1 


< 


1+q 2 _ 



1/2 



Thus for all values of q in 
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Proof: To begin with we write 



4>l+i = ^Pt 1 ® IT) + v^-Vi ® |0) 



(2.4.81) 



where is a spin-1 groundstate on a chain of length N with total S 3 component 
M. Define 



-L-l 



IFl II 



Now ex = ^—^prr^-) so we calculate G[l,l+i\^> 



G 



[L,L+l]</> 



G 



^^^- 1 ®IT)-^®|0) 



Ittt^I? ® |m = 2) + 



II^M 2 



1 



2g 



-21 



®\m = 2) 



L -!||2 



- - 1 Vr i ® |m = 1) 

v/TT^Vll^ 1 !! 2 



L-2M2 

m^ L L-l ®\m=l) 



where |m = k) is the ground state of i?[L,L+i] such that 

(££ + ££ + i)|m = fc) = fc. |m = k) 
Calculating ||G [LiL+ i]9?|| 2 , we get 



l|G[L ; 



L+l 



II^M 4 



Ut-l\ \ A 

l + q 2 



2q 



- 2L - 2 u L L zi\\ 2 + 



But 2q~ 2L ~ 2 + 2q 



-2L 



On the other hand 



H^Iill 2 (2g- 2i " 2 + 2g" 2i 



llv>M 4 



«ll 2 



l + g 2 
I 2 so we have 



\\G[L,L+lM\' 



«H 2 | 



^ + ll| 2 



ii^r i ii 4 (i+? 2 ) 



;£ + J 2 



^11 2 



«ll 2 ) 



(2.4.82) 
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so that 

4= llG[^ll a = ^g. 1 (2 . 4 . 83) 

IMP ||^M 2 i + <? 2 1 ; 

Substituting ||Vi _1 || 2 = 2^ =1 f 21 into (gXHj gives the expression 

, (1 ~ q 2L ) q 2 , , 

£l ~ (i - ' IT?' (2A84) 

Hence e| < 1/2 since for g e (0,1), (1 - g 2L )/(l - g 2L+2 ) < 1 and g 2 /(l + q 2 ) < 1/2. 



What remains in order to verify assumption (|2.4.fiHj) . is to show that in sectors 
with m < L, e m < e^. As stated previously this is really the hard part. At this point 
we are unable to prove this analytically, but we state the conjecture here and provide 
a little justification, both numeric and heuristic. 

Conjecture 2.4.8 e m < el for all m such that —L < m < L. 

Some of the strongest evidence we have for this conjecture is numeric. We show 
here the values of e m in all sectors for a chain of length eight and various values of A. 

The general shape persists for all values of L that we have tried (up to 11,) and 
moreover, for fixed A the values of e m stabilize very quickly relative to the number 
of downspins, i.e. in the sector corresponding to k downspins for chains of length six, 
seven, or eight the SL-k varies no more that one part in a thousand. We include our 
code at the end of this paper in Appendix O 

It is easy to compute that in the Ising limit (A — > oo), we get e m = for all 
m. We hope to prove Conjecture 12.4.81 for at least most values of A by performing 
perturbation around the Ising case. This should be possible since for long chains the 
ground states of the XXZ Hamiltonian can be thought of as perturbations of the Ising 
ground states. One also notices that e m becomes very small, i.e. close to the value 
in the Ising case, for relatively small values of A (A « 2). This behavior gives us a 
hope of being able to prove Conjecture 12.4.81 for at least some values of A. 
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Figure 2.1: Some plots of e m for a chain of length L = 8 for different values of A. 

The following is a plot comparing the true gap of the XXZ Hamiltonian for a chain 
of length L to the lower bound given by Theorem 12.4.51 The anisotropy parameter 
A varies between 1 and 20. 

The plot shows that the estimate captures the behavior of the true gap very well. 
Over the inteval 1 < A < 20, the estimate is only off by 17% on average. The 
maximum error occurs for A close to 1 and at this point the estimate is off by about 
70% from the true gap, but as A grows and gets closer to 20 the error shrinks to 
about 7% of the actual gap. 
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A comparison of the true gap ot the XXZ Hamiltonian vs. the 
lower bound given the Martingale Method. 
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Figure 2.2: The upper plot is the true spectral gap of the XXZ chain on L sites. The lower plot is 
the estimated gap from Theorem 12. 4. 51 (1 < A < 20) 
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Chapter 3 

XY Model 

3.1 Introduction 

Understanding phase transitions, for example ice to water, or non-magnetic to mag- 
netic, is one of the central questions of statistical mechanics. Spin systems were 
initially introduced by Lenz in 1920 to try to explain the tendency of some materials 
to possess macroscopic magnetic fields. He proposed a simple model in which atoms 
in a crystal could assume two possible orientations which later came to be called 
the Ising model jSj. It was Ising's failure, in 1925, to discover a phase transition in 
Lenz's model which led Heisenberg to introduce quantum spin models in 1928. In 
1935, Peierls gave an argument that demonstrated that the Ising model does exhibit 
phase transitions in two and greater dimensions [21] • Although Peierls' argument was 
not completely rigorous it did contain all the essential ideas required for the desired 
result. The Peierls' argument was later made rigorous by Griffiths jU], and Dobrushin 

Mathematically, phase transitions are associated with a singularity in some ther- 
modynamic potential, generally the free energy. While the Ising model is a classical 
model, it has become the basis for understanding magnetic phase transitions in quan- 
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turn mechanical models. The Ising Hamiltonian is given by 



h a({ s x}) = ~ ^2 s * s y- h ^2 s x (3.1.1) 

x,y£A XgA 

\x-y\=l 

where A is a finite subset of Z d , h is some positive number, and {s x } is in the set 
Q\ = {±1} A . The partition function of the system Zj^(h, (3) is given by 

z{{h,p)= Yl H i(M) (3.1-2) 

{s x }en A 

For any thermodyamic variable, X, defined on Q\, we define the expectation of X at 
inverse temperature (3 by 

In dimension two and greater the Ising Model exhibits spontaeous magnetizaton. 
Consider the average magnetization at a particular site, m(h, (3) = uJh,/3(so). Then 
there is a constant [3 C = /3 c (d), such that for all j3 > j3 c , m(h, (3) is discontinuous at 
h = 0. Specifically 

lim m(h,(3) = - lim m(h,/3) (3.1.4) 

for all (3 > (3 c (d). This implies that the free energy, f J (h,f3), is not analytic at h = 
since 

m(h J (3) = ^j-(f I (h J (3))\^ (3.1.5) 

Physically, the discontinuity in m(h, j3) at h = 0, can be interpreted as the system "re- 
membering" the effects of an external magnetic field after the field has been removed, 
i.e. it continues to possess a macroscopic magnetic field. 

The phase diagram of the XY model has been the subject of much study. For high 
temperature, one can use a high-temperature expansion to show that there is a unique 
phase. To see an excellent discussion of such techniques see j2Ej- In [Zj Dyson, Lieb, 
and Simon showed that, in the absence of an external magnetic field, the spin-1/2 
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XY model does have a phase transition for positive temperature in dimensions greater 
than two. The Mermin- Wagner |3] theorem shows that for positive temperature this 
result can not be extended to dimension two. In a later paper ^3], Kennedy, Lieb, 
and Shastry demonstrated that in dimensions two and greater the ground state of 
the XY model has off-diagonal long-range order. 

T = 1/p 



1 




Figure 3.1: The phase diagram for the XY Model. The numbered regions are: 1. (h = 0) For 
dimension three and greater Dyson, Lieb, and Simon showed that LRO exists for T > 0. 2. For 
spins bigger than one-half and dimensions two or greater, Kennedy, Lieb, and Shastry showed the 
groundstate has LRO. 3. The low-temperature/small h regime is the most interesting, but nothing 
is known (the shape of the curve is unknown.) 4. We show that there is a low-temperature/large h 
region in which there is no phase transition. 

Here we show that the free energy of the XY model in the presence of an exter- 
nal magnetic field above a critical strength, is analytic for some region of non-zero 
temperature in dimension two and greater. The more interesting regime is when 
h < h c = 2d, where there is a connection to Bose-Einstein condensation [T^]. Our 
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results do not extend to this case. 

We consider a perturbation of the spin-- XY model by a uniform magnetic field 
in the z direction. The strength of the magnetic field is great enough so that there 
is a unique ground state of all "up" spins. We use a contour expansion inspired by 
Kennedy in [X2| to extend the uniqueness of phase to non-zero temperatures. 

/ \®|A| 

Take 7i A — ( C 1 where A C Z , d > 2, and |A| < oo. For x, y e A such that 
\x — y\ = 1 we write {x, y} = (x,y). The Hamiltonian is defined as follows: 

H' A = H^ Y + hHf, h>h c (d) = 2d (3.1.6) 

where 

(x,y)cA (x,y)cA 

and 



The operators a 1 , a 2 , and <r 3 are the Pauli spin matrices 





a" 



1 

10/ \ i I \ -1 

where the subscript denotes the site on which they act. We take the strength of the 
magnetic field h to be greater that some h c (d) so that the unique state is the all up 
state ||) A given by 

ir> A = (gm>, (3.i.9) 

For the remainder of this chapter we will work with the scaled Hamiltonian 

H A = (l-5)H XY + h c H z (3.1.10) 
where the parameter 5 is given by the equation 
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By (j3.1.6j) . h > h c so we have 5 G (0, 1). 



Let 



Z A (<*,/?) = Tr e~ pHA 



(3.1.12) 



be the partition function. Our goal is to show that the free energy 



/(*,£) = -! lim ±log[Tre-^ A ] 

p A/"Z d \A\ L J 



(3.1.13) 



is analytic for all (5, (3) in some appropriate region, and hence that no first-order 
phase transition exists. In order to do this, we will express Z\{6, 0) as a contour 
expansion. Then, we follow Kennedy in [12], and show that the partition function 
can be rewritten as a cluster expansion. The key to this is Theorem 13.3.21 which 
shows that the weights of the contours decay in an appropriate way. Once the cluster 
expansion is established the analtyticity of f(5,/3) follows immediately, as the cluster 
expansion is a uniform limit of analytic functions. 

3.2 Basic Definitions 

A set A C A will be called a contour. This differs from the usual definition of 
countours as set of bonds separating lattice sites with opposite spins. This is okay 
though, because in the model we consider the ground state does not possess the spin- 
flip symmetry that generally introduces ambiguities. A connected set P C A will be 
called a polymer. To each contour A G A let \ A) denote the basis vector of H,\ such 
that 



Furthermore, if (o"+cr +a x a+)\A) ^ then we will denote the corresponding contour 




1.4)* = |T) xgA\^1 
14)* =|l> x^A 



(3.2.14) 



by + o x a 



)A. 
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Returning to Z A (5,j3), by the Trotter product formula, we can express e ^ H as 

N 



e -l3H A = j im 
N — >oo 



0h c ttZ I 3(1 — 5) ttXV 



so that 



Z A (6,P) = Tr e~ pHA = lim Tr 

N — >oo 



_£hc H Z 

6 N H A I i 



P0- - *) U XY 



N 



m 



N 



hm e * A ( 1 

JV >oo ^ V ' V 

aCA L v 



iV 



if 



2V 



\a) 



inserting a resolution of the identity between each term in the product we can write 
the above as 

8(1-5) 



Xre -/3H A = lim y- / 

N — >oo t-^ 



«0 



• • • (ttjV-l 



— H A 1 



AT 



Z^ 1 - *)jyxy 



e n h a i _ 



iV 



-if 



K) 



XY 



|ai)(«i| • • • 



(3.2.15) 



Taking this new expression for Z\(5, 0) we see that the non-zero terms in the sum are 
characterized by sequences of contours {a , a%, . . . , «tv} with the following properties 

(i) a = a N 

(ii) a i+1 = < or 

(o"+CT~ + a x a y) a i f° r some ( x , U) 

This motivates the following definition 



Definition 3.2.1 A quantum contour Y = {r ,r 2 , 
tours satisfying 



T N } is a sequence of con- 



(i) r = r 



A 1 



(ii) r n+ i = < 



or 



{a+a + a x a+)r n where \x — y\ — 1 
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Figure 3.2: The spin-flips for the above quantum contour are 1 — ► 2 : (CA, DA), 2 — > 3 : (FA, GA), 3 
4 : (F3, G3), 4^5: (F3, FA), 5 -► 6 : (G3, GA), 6 1 : [CA, DA) 



We now write 



Za{6,P) 



lim 

N — »oo 



ao,...,o/v-i 



e ~— H A 1 



0(1 - %XY 



N 



lim Y^cufr) 



N 



K) 



where 



w(7) 



-ph c 



r„| W~*) 



n(T) 



|ai)(ai| 



(3.2.16) 



(3.2.17) 



N 

and n(T) is the number of times r i+ i = {p^a~ + <7~er+)rj in the quantum contour T. 

Let S(T) = |jr n be called the support of a quantum contour T. If S(T) is 
connected, then T will be called connected. We will use 7 to denote a connected 
quantum contour. Two connected quantum contours 7 1 and 7 2 are disjoint if S^ 1 ) U 
S'(7 2 ) is not connected. 
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+++++++++ 



Figure 3.3: The support of the quantum contour above. Note that it is connected. 

For a given quantum contour T, S(T) can be written as the disjoint union of 
connected parts 

k 

S(X) = \_\S*(T) 
i=i 

Then, T can be split up into connected parts by 

T £ = r n n^(r) 

and we write T = Uryi . We call two general quantum contours r\ and T 2 disjoint if all 
there connected components are disjoint in the sense of connected quantum contours 
defined above. 

With these definitions fixed consider a quantum contour T = T 1 U T 2 , then a>(T) 
will factorize in the following way 

n(T) 



But, T = T 1 U T 2 so that 



w(n=e w«izi>y 



|r | = K\ + \Tl\ (3.2.18) 
n(T) = n^ 1 ) + n(r 2 ) (3.2.19) 



Thus 



u(T) = e-^^^j' e-^ol^^y) (3.2.20) 
^(r 1 )^^ 2 ) (3.2.21) 
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So that 

Tre -/3H A = lim y^uCr) 

N — >oo ^— ' V ' 

r 

oo 1 n 

£ (3.2.22) 



TV — >oo z — ' n! 

n=0 ~1 ~n j = l 

all disjoint 



Now, for a polymer P, we define 



W(P)=lim W W ( 3 - 2 - 23 ) 

iv — >oo ' 

7 :<?( 7 )=P 

Using these new weights we can re-write the partition function as 

oo 1 n 

Z A (<5,/3) = Km £^ £ n<"M ( 3 - 2 - 24 ) 

n=0 -.1 . . ~n J = l 

all disjoint 

oo n 
n=l {P 1 ,...,P„} i=l 

These are all the basic definitions which we need to apply the cluster expansion 
3.3 Cluster Expansion 

Cluster expansions were introduced as a perturbative method for high-temperature 
gases. They apply to continuous systems, e.g. classical interacting particle systems, 
as well as discrete systems such as lattice particle models, polymer models, or, as in 
our case, spin models. There are many good references, perhaps the most accessible 
being Simon's book on lattice gases [21], while the treatise by Brydges jlj is also 
excellent. Others include [101, [IB!) E3' EL 0, and [20] The basic goal of the cluster 
expansion is to give an expression for the partition function of some system that can 
be more easily manipulated that the standard expession as the sum of exponentials. 
Here we present a condition for the convergence of the cluster expansion attributed 
to Kotecky and Preiss [TH]. The proof is taken from Ueltschi in [3J. Although the 
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setting of (3T] is both continuous and discrete space, here we only need the discrete 
case. 

Let A be a finite set whose elements are called polymers. Let t be a symmetric 
and reflexive relation on A x A. For polymers A,A'eA we say that A and A' are 
incompatible if AiA' and are compatible otherwise. The partition function of the 
polymer model is 

oo n 
n=l {A lt ...,A„} i=l 

where the to is a complex function on A and the sum is over sets of polymers whose 
elements are pairwise compatible. 

A sequence C = (Ai, . . . , A n ) of polymers is called a cluster if {A\, . . . , A n } can 
not be partitioned into two sets, all the elements in the first being compatible with 
all the elements in the second. We can make this definition more visual. Consider the 
graph G(C) with n vertices where there is an edge between % and j whenever A{iAy 
C is a cluster iff Q(C) is connected. We write G C G{C) if 6? is a graph with the same 
number of vertices as G{C) and whose edges form a subset of the edges of G(C); \G\ 
denotes the number of edges of G. Define a function, <p(C), on clusters given by 

f 1 if n = 1 

<P(C) = < (3.3.27) 

[ X]GcS(C*),conncctcd( — 1) if ^ > 2 

and a function $(C) on clusters by 

1 n 

$(C) = - V (C)Hw(A t ) (3.3.28) 

n! ±J - 

i=i 

where C = (A\, . . . ,A n ). The following theorem is due to Kotecky and Preiss |18j . 
It says that if the weight function on the polymers, w, satisfies some condition then 
the logarithm of the partition function Z(A) can be given by a sum over clusters. 

Theorem 3.3.1 Let a and b be non-negative functions on A such that for all A £ A 

\w{A')\e< A '^ h W < a{A) (3.3.29) 

A'A'iA 
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Then, 



Z(A) = exp HC)j (3.3.30) 



where the sum converges absolutely. Furthermore 

J2MC)\e b{C) < a(A) (3.3.31) 

CiA 

where CiA means that A is incompatible with at least one polymer inC = (A\, . . . , A n ) 
and b{C) = YJU KA). 

So, in order to show that f(S,(3) is analytic, we prove that the weights of the 
polymers W(P) given in ()3.2.23|) satisfies condition (|3.3.29j) . Then Theorem 13.3.11 
allows us to write 

f(6,/3) = ~ lim ±rhxZ A (6,l3) (3.3.32) 
p K/z d | TV | 

= -- lim i- V$(C) (3.3.33) 
P k/^ |A| ^ 

which is the sum of analytic functions. Then the analtycity of f(8,P) follows from 
the uniform convergence of this sum. To this end we state and prove the following 
proposition and leave the proof of Theorem 13.3.11 to Appendix El 

Proposition 3.3.2 Given 5 where < S < 1, there exists a Pq and an e such that 
for all (3>(3 

J2w{P)e £lPl < e (3.3.34) 

Ply 

where y G Z 2 is fixed. 

Propostion ()3.3.2j) immediately gives us ()3.3.29j) since 

W{P)e £ \ p \ <J2J2 W{P)e £ \ p \ <J2 £ - £ \ P '\ (3.3.35) 
p-PnP'^% yeP' P3y yeP' 

Proof: Initially we will express the sum Y^P3y W(P)e £ \ p \ as a sum over contours. 
Recall that 

W(P)=hm "(7) 
r s(ri)=P 
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So, J2pb v W(P)e e ' p ' can be re-written as 



Y^w(p)e^ p \ = £ him T, w W eE|p| ( 3 - 3 - 36 ) 

Pay P^y \ r -S(r/)=P J 

= lim V cu( 7 )e £|5(7)l (3.3.37) 

where we have written 7 : 7/ for the condition y G ^(7). Let the new Hamiltonian 
be given by 

H° A = (1- 5/2) (H* Y + h c Hl) (3.3.38) 
and define the modified weights ujo(j) by 

n(7) 



0,0(7) = e -«i-wd7oi Ai-V2)/^ y 



(3.3.39) 



which are just the weights of a quantum contour 7 when i?A is replaced by H\ in the 
expansion of the partition function Z\{8,(3) The weights u(j) and ^0(7) are related 
by 

wiri) = m(i)uo(y) (3.3.40) 

where the multiplyer 772(7) is 

m ( 7 ) = e -/'fAclTD| p »(7) j p = 1 ~ * < 1 (3.3.41) 

1 — 0/2 

Thus, the sum on the right side of equation (|3.3.37|) can be expressed as 

£ ^( 7 )e £|5WI = £ ^ (7)e- /3 ^ c|7ol p n(7) e £ l 5 WI (3.3.42) 

7 :j/ 7:5/ 

Let -^(7) be the set of spins that flip at least once in 7 then ^(7) satisfies the 
inequalities 

\S(^)\ < |F( 7 )| + | 7 o| (3.3.43) 

and (3.3.44) 

|F( 7 )| < 2n( 7 ) (3.3.45) 
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Inequality (J3.3.43|) follows since ^(7) contains all the spins that flip and any 73 con- 
tains all the spins that never flip, and we get inequality (|3.3.45|) from the fact that 
for every occurence of a spin flip in 7 two spins change, but it might be possible for 
certain spins to change more than once in the course of a quantum contour 7. 

From inequalities ()3.3.43j) and (J3.3.45|) we get |5'(7)| < 271(7) + ] To i , which in turn 
implies 

e ^\S(j)\ < e 2en(7) e e|7o 

From the above inequality we can determine 

for (3 large enough, and e small enough. Applying this estimate to ()3.3.42|) we get 

oj{l)e^ < ^o(7)p^e-^ h ^°l (3.3.47) 
r-y r-y 

We now re-write the sum on the right side as a sum over initial configurations 70 and 
then a sum over quantum contours with that intial state. This give the expression 

£ uMp^e-*^ = E E ^o(7)p^e-^l G l (3.3.48) 
r-y GcAr-G,y 

where we have used 7 : G to indicate the sum is over quantum contours 7 such that 

7o = G. 

Let D y (G) be the number defined as follows: A set K C A will be called G,y- 
connecting if the set (K U G U {y}) C A is connected. D y (G) is then defined as 

D y {G) = min{|J-s:| : K is G,y - connecting} (3.3.49) 

In other words, D y (G) is the minimum number of sites required to make G U {y} a 
connected set. Since 7 is a connected contour, ^(7) is 70, y-connecting, hence 

D y ( T o)<|F( T )|<2n( 7 ) (3.3.50) 
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This allows one to bound the term p n ^l 2 by p D y\ia)/ A . The advantage to this is that 
now the term does not depend on the choice of quantum contour, but only on its 
initial configuration. Therefore, in ()3.3.48|) . we can write 

£ £ c ( 7 )p^e-^ic| < £ p ^ e -^\c\ £ , o(7) (3 .3. 51) 

GcAr-G,y GcA 7:G 

We can drop the condition that 7 be connected since adding more terms only increases 
the sum. Thus, by ()3.3.47j) . (l3.3.48j) . and (|3.3.5lj) . we have 

lim Vcj( 7 )e £ l 5 WI< lira r^e-^ |G| Vw (r) (3.3.52) 

N—>oo ' J N-^oo ' J ' J 

r-y GcA r-.G 

However, the only term on the right side of (j3.3.52j) which depends on the length of 
the quantum contour T is 

T:G 

thus we can pass the limit through the initial sum to get 



lim V uj {T) = Tr(e- pH ^o¥ G ) (3.3.53) 
r-.G 

= e-^<^ (3.3.54) 

where P G = \G)(G\. But H° A > 0, so e - ^^ < 1 Thus we have shown 

lim Y a;( 7 )e £ l^>l < V p^e^t^l lim V cu (T) 
r-y GcA r-.G 

< y^p^e-p^m (33 55) 



f 

GcA 

The proof is then completed by the following lemma 



Lemma 3.3.3 For fixed 5, < 5 < 1 

Dy(G) „ S 

f 

GcA 

where r(S, (3) — ► as (3 — ► 00 



yp^e-^Ur&P) (3.3.56) 
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To prove Lemma 13.3.31 we follow Kennedy and think of G as being generated 
by a random walk on points on the lattice in the following way: Let F be a G, y- 
connecting set with \F\ = D y (G). Consider a nearest-neighbor walk fj on the elements 
of F U G U {y} where fjo = y- We take f/ in such a way that it does not visit any 
site in F U G U {y} more cq times, where Co is a constant that only depends on the 
dimension. We take G to be generated by a walk rj where 77 is the walk f] with the 
steps on the set F removed. Notice that in general 77 is not a nearest-neighbor walk. 
The product 



Dy(G) „ s 



p 4 e 



(3.3.57) 



contains a factor p 1 / 4 for every site in F and a factor e ^ c for each site in G. Taking 



gives 



Po = P 4c ° an d S 



Sh c 

4Cn 



p 4 e r < 



Po 



(3.3.58) 



(3.3.59) 



i=l 



where d(r]i_i,r]i) is the minimal number of sites required to connect r]i_i and 7/j. 
Therefore 



But 



Po 



GcA 



'/ 

00 



8=1 



Ell"! 



n=l 



r):\rj\=n i— 1 



(3.3.60) 
(3.3.61) 



?7:|rj|=n *=1 



d(»7i-l > »?i) 



e n« 

6 ,...,6 n eA i=0 



d(bi-i,bi 



< 



E" 



d(6 ,6) 



(3.3.62) 
(3.3.63) 
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So, let K = supv J2bPo^ b ■ Then 

oo 

p^e-^^' < K n e~^ n = r(S, (3) (3.3.64) 

GcA n=l 

and we note that r(5, (3) — > as (3 — > oo. 

This completes the proof of Lemma 13.3.31 which in turn completes the proof of 
Propostion 13.3.21 
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Appendix A 

Calculation of critical magnetic 
field. 

In Chapter|3]above we presented a contour expansion for the free energy of the system. 
In this expansion the contours that contribute the most are those that are close to the 
ground state of the Hamiltonian. They can be thought of as islands of down spins in 
a sea of up spins. The argument relies upon the fact that we only consider a region in 
which the all-up spin state is the unique ground state of the system. Here we present 
the proof that there is a critical magnetic field strength, h c = 2d, such that for all 
h > h c the all spin up state is the unique ground state of the XY Hamiltonian 

H XY = - E + O + h £ \& ~ 

\x—y\=y x 

We show h c > 2d by calculating the energy of spin waves and comparing it to the 
energy of of the all-up state. The inequality h c < 2d is obtained by writing the 
Hamiltonian as the sum over two-site interactions and then determining a condition 
on h such that the all up state minimizes every two-site interaction. 

Proposition A. 0.4 There is a critical value h c such that for h > h c the ground state 
of Hi is the all 'up' state, while for h < h c this is no longer the case. Moreover 
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h c = 2d 

Proof: To find the critical value of h we first calculate when a spin wave has lower 
energy than the all up state. Let ip k = Y2 X a x&~ -|t) where a x = e^'/IA] 1 / 2 . Then 



{^ k \Htij) k ) = (%l) k \H xx %l) k ) + (tfj k \H z ij k } 

-a + )a-,m + 



z,z' z > z 



But (tjcr+cr+cr y cr^lf) = 5 x ^S Vt g and (11°^°* a y a z ' It) = ^.A,*') so 

{il) k \Hi) k ) = -^a y a x -^a x a y + h (A.0.1) 

= -2^^ + /! (A.0.2) 

2 1 d 

= E e ~ lkX 2 J2(e ik{x+ej) + e ik{x - e > ] ) + h (A.0.3) 
' ' xeA j=i 

d 

= -2^cosA; i + h > -2d + h (A.0.4) 

3=1 

Therefore the spin wave ip k at least has variational energy — 2d + h. But the state ||) 
has energy zero, so we are looking for the regime in which 

-2d + h<0 ^ h<2d (A.0.5) 

This implies 

h c > 2d (A.0.6) 

To get the other inequality we write 

H = E (#2 + Yd^j ' ^-2 = ~ {a> y + a * a y ] ~ Td {al + a 'y ] (A -°- 7) 

(2) 

If we denote by E and e the ground state energies of H and respectively, then 



we see 



£o>5>° + |p (A.0.8) 
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(2) 

Then we can calculate for which values of h the state |f) minimizes each Hx,y. For 
these values of h 

E = £(e„ + A) 

The eigenvalues of H^l are ±1 and ±h/2d. The eigenvalue —h/2d corresponds to 
the state Iff) a t the sites x and y. Hence, for ||) to minimize H^l h must satisfy 

-— < -1 /i > 2d (A.0.9) 

which gives the inequality 

/i c < 2d (A.0.10) 
Therefore, we have 2d < h c < 2d which implies h c = 2d. | 
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Appendix B 



Proof of Theorem 13.3.1 



We restate Theorem 13.3. II here and present a proof due to Ueltschi The proof is 
very technical and is included here for the sake of completeness. Let A be a finite set 
whose elements are called polymers. Let t be a symmetric and reflexive relation on 
A x A. For polymers A, A' 6 A we say that A and A' are incompatible if AiA 1 and 
are compatible otherwise. The partition function of the polymer model is 



where w is a complex function on A and the sum is over sets of polymers whose 
elements are pairwise compatible. 

Theorem B.0.5 Let a and b be non-negative functions on A such that for all A £ A 



oo 



11 



n=X {Ai,...,A n } i=l 




(B.0.1) 



A',A'lA 



Then 




(B.0.2) 



where the sum converges absolutely. Furthermore 




(B.0.3) 
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where CiA means that A is incompatible with at least one polymer inC = (Ai, . . . , A n ) 
and &(C) = £r=i b{Ai). 

The sequence of the proof is out of order with the statements of the theorem. We 
first prove inequality (jB.0.3|) as this follows from the assumption (jB.0.1|) . We then 
use flB.0.3|) to prove the convergence of the cluster expansion on the right hand side 
of ()B.0.2|) and finally show that the partition function Z(A) can indeed be written as 
the cluster expansion. 

Proof: In order to prove inequality (|B.0.3|) we induct on number of elements in the 
cluster. For n — 1, C — (Ai) then by assumption (|B.0.1|) we have 

MC)\e b( ~ c) < a(C) 

ClA,\C\=X 

Assume that for all k < n 

|$(C)|e fe(c) < a(C) (B.0.4) 

dA,\C\<k 

and consider the sum over clusters with n + 1 elements or less. To bound this, first 
sum over polymers A\ incompatible with A and then over remaining polymers. Thus, 
we have 

n+1 1 

ClA,\C\<u+1 A x lA j=2 U > A 2 ,...,Aj 



l\w{At)\e b ^\ Yl (" 1 )' G| ( B '°- 5 ) 



Gcg(A 1 ,...,A j ) 
connected 



For a given graph G, let (Gi, . . . , Gk) be a sequence of connected subgraphs of G. Let 
each subgraph Gi have vertice set V{ where V^H Vj = if % ^ j. Furthermore, assume 
that V\ U • • • U T4 = {2, . . . , n}. The sequence (G±, . . . , G n ) defines a disconnected 



APPENDIX B. 



Convergence of Cluster Expansion 



58 



subgraph G' of G obtained by removing all edges eminating from vertex 1. Hence, 

E (-^'I^EfI E ni-^D-i)™ 

acg(A 1 ,...,A j ) k>l (Gi,...,G fe ) i=l G^ 

connected 

=EfI e n(-i) |o "| (B.o.6) 

fc>l (Gi,...,G fc ) i=l 

since Y^ G / (— 1)' G * = —1. Next we expand the sum over sequnces (6?i, . . . , G^). First 

i 

sum over the respective number of vertices mi, . . . , with the condition that m\ + 
■ ■ ■ + rrik = j — 1, and then choose connected graphs for each choice of mi, . . . , vrtk- 
We can bound the sum over A2, ■ ■ ■ , Aj by summing over clusters of length k that are 
incompatible with the polymer Ai such that |Cj| — mi. This gives 

E l*(C)|e" C| <E»^)e^5:Es 

A- 



CiA,\C\<n+l 



A 1 lA 



j=2 k>0 



e e ni*wi 

(mi,...,m fc )hj-l Ci,...,^ j=1 
GitAi,|Ci|=m,i 



(B.0.7) 



Relaxing the condition that (mi, . . . , m^) \- j — 1 < n to rrii < n for all i and then 
summing over j gives 



MC)\e bic) 

Ct,A,\C\<n 



(B.0.8) 



-1 fe 

Cti4,|C|<n+l AiiA k>0 

We now use the induction hypothesis (|B.0.4|) to bound the term in the brackets by 
a(Ai), which in turn bounds the sum over k by e a ^ Al \ But, under the assumptions 
of the theorem (jB.0.1)) 



Y \w(A 1 )\e b{M)+a{M) <a(A) 



(B.0.9) 



A11A 



which give (jB~fT8|l by induct ion. 
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To get the convergence of the cluster expansion in ()B.0.2|) note that any cluster 
is incompatible with its first element, so 



Ei^i^EE^^ 00 



(B.0.10) 

C AeA ClA 

by (jB.0.3|) . Thus we have established the convergence of the cluster expansion, what 
remains is to show 

Z(A) = exp (j2 $(C) 

This will be done by expanding Z(A) and correctly rearranging terms so that the 
sum over clusters becomes apparent. 



Z(A) 



E U w ^ 

compatible 

1 n 

1 + E^ E n^IK 1 -^^]) 



n>l (Ai,...,A n ) * =i 



! + E^ E n«(*) E 



-1) 



|G| 



(B.0.11) 

(B.0.12) 
(B.0.13) 



n >! (Ai,...,j4.») i=l 



GC6(A 1 ,...,A„) 
connected 



Again we sum over partitions of the vertices m = (m 1; . . . , rrik) \~ n, first summing over 
the size of the partition k, considering sequences of graphs (6?i, . . . , Gk) corresponding 
to the partition m. This gives 

^) = i + EEA E — ^ — r 

n>l fc>l (mi,...,mfe)hn 



n 

i=l 



e n«w e ( _i ' igi 

J 4l,...,Am i J = l 



Gc5(A 1 ,...,A mi ) 
connected 



(B.0.14) 



1+ EE^ E II E *w 

n>l fc>l (mi,...,mfe)hn i=l C,|C|=mj 



exp ^ $(C) 
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where the last line follows since by (jB.0.10|) we have absolute convergence of the 
cluster expansion. g 
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Appendix C 

Matlab Code 

We include the Matlab code with which we produced the images in Chapter |3 There 
are two main programs spinlkink.m and spinlcompare.m. In spinlkink.m we evalu- 
ate ||G[l,l+i]-Ex|| for L+ 1 = 8 in each sector of Sf ot . This is accomplished by building 
the XXZ Hamiltonian using the spin operators Sz, Sx, and Sy defined below. The 
Matlab function kron computes the tensor product of its inputs. The projections 
onto ground states are produced by calculating ground state vectors and then creat- 
ing the operator \v)(v\ where v is the given ground state. Finally we compute the 
proper eigenvalues and plot them versus the Sf ot sector. 

spinlkink.m 

%spinlkink.m 

% SPIN 1 KINK HAMILTONIAN 
clear; 

delta =50; q = delta -sqrt (delta~2-l) ; 
A = sqrt(l-delta~(-2)) ; 
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G2 



Sz=[l 0;0 0;0 -1] ; 

Sx=[0 1 0;1 1;0 1 0]/sqrt(2); 

Sy=[0 1 0;-l 1;0 -1 0] /(i*sqrt(2)) ; 

hdiag= kron(Sz,Sz); 

hhop = kron(Sx,Sx) + kron(Sy,Sy) ; 

hbdry= kron(Sz,eye(3)) - kron(eye(3) ,Sz) ; 

HDIAG = sparse(3~L,3~L) ; 

HHOP = sparse (3~L, 3~L) ; 

HBDRY = sparse(3~L,3~L) ; 

Sztot = sparse(3~L,3~L) ; 

for a=l: (L-l) 

HDIAG = HDIAG + kron(speye (3" (a-1) ) , kron(hdiag, speye (3~ (L-l-a) ) ) ); 
HHOP = HHOP + kron(speye(3~(a-l)) , kron(hhop, speye (3~ (L-l-a) ))) ; 
HBDRY = HBDRY + kron(speye (3" (a-1) ) , kron(hbdry, speye (3" (L-l-a) ))) ; 
Sztot = Sztot + kron(speye(3~ (a-1) ) ,kron(Sz, speye (3" (L-a) ) ) ) ; 
end 

Sztot = Sztot + kron(speye (3~ (L-l) ) ,Sz) ; 

Hkink = - HDIAG - 1/delta * HHOP - A*HBDRY; 
H = Hkink + (L-l) * speye (3~L) ; 

°/„ Projection G onto the last two sites, [L,L+1] 
h2 = -hdiag-l/delta*hhop - A*hbdry + speye (3~2) ; 
[v,d] = eigs(h2,5, ' sm' ) ; 
G=0; 

for k=l:5 



APPENDIX C. Matlab Code 



G = G + v( : ,k) *transpose(v( : ,k) ) ; 

end 

G = kron(speye(3~(L-l)) , G) ; 

for downspins = 1:2*L-1 
for m = downspins-1: downspins+1 
Proj =[] ;Proj = speye(3~L); 
for n=0: (2*L) 
if ne(n,m) , 

Proj = Proj*(Sztot - (L-n)*speye(3~L))/(n-m) ; 
end; 
end; 

Hred = transpose (Proj )* (H - 0.2*speye(3~L))*Proj 
[V,D] = eigs(Hred + . 2*speye(3~L) , 1 , ' sm' ) ; 
psi ( : ,m-downspins+2) =V; 
end 



Psi(:,l) = kron(psi(:,l), [0,0,1]'); 
Psi(:,2) = kron(psi(:,2),[0,l,0]'); 
Psi(:,3) = kron(psi(:,3),[l,0,0]'); 



A=[] ; 
for i=l:3 

for j=l:3 

A(i, j)=transpose(Psi( : ,i))*G*Psi( : , j) ; 

end 

end 
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G4 



epsi = eigs(A); e( : ,downspins+l) = epsi(2); 
end 

e(:,l) = sqrt(q~2/(l+q~2) * (l-q~ (2*L) ) / (l-q~ (2*L+2) ) ) ; 
for i=l:2*L-l 

p(i)=e(2*L-i) ; 
end 

s=[-L+l:L-l] ; 
plot(s,p, 'b:*') ; 

The program spinlcompare.m is very straightforward. It calls two other programs, 
Hmltn.m and SectrProj.m. Hmltn.m produces the XXZ Hamiltonian of inputted spin 
on a chain whose length is passed to the program. SectrProj.m produces a projec- 
tion onto the desired sector of S^ ot . In spinlcompare.m we loop through values of A 
between 1 and a specified number DeltaMax. For each value of A spinlcompare.m 
calculate the spectral gap of the output from Hmltn.m by projecting into each sector 
of S" t 3 ot and calculating the two smallest eigenvalues. In each sector, is one of the 
eigenvalues and the other is the spectral gap. spinlcompare.m then plots the spectral 
gap for each A along with the estimate from Theorem 12.4.51 

spinlcompare.m 

% spinlcompare.m 

°/ This code produces a graph that compares the true gap of the 
%XXZ Hamiltonian with the estimate in this Thesis: Theorem 2.4.5 
°/ The variables are: 

% L = Length of the chain J = spin of the system considered 

% DeltaMax = this is the max value of Delta (anisotropy parameter) 
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clear; 
L = 8; 

J = i; 

DeltaMax = 20; 



for i=10: (10*DeltaMax) 
Delta = (i/10) ; 
H = Hmltn(L, J, Delta) ; 



°/„Sets H to be the XXZ Hamiltonian 
°/ on spin-J chain of length L with 
°/>anisotropy Delta 



for j = -L+1:L-1 

P = SectrProj (L, J, j) 



°/ This sequence calculates the 
7»gap of H by projecting into 
/^different sectors and then 
smallH = transpose (P)*H*P; °/ finding the two smallest e-value: 

°/„Since there is a unique 
"/.ground state in every sector the 
%first non-zero eigenvalue 
°/ is the gap in that sector. 



v = eigs(smallH,2, 'sr') ; 



u(L+j) = v(2); 
smallH= [] ; P = [] 

end 

gap(i-9) = min(u) ; 



%This asks now for the minimum 
°/ gap over all sectors 
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gamma = 2.5 - sqrt(2.25-(2/Delta~2)) ; %This is the gap of the 

°/>two-site Hamiltonian gi 
°/ in 2.4.75 

q = Delta - sqrt (Delta~2-1) ; °/ q solves equation 

7 2Delta = q + q"(-l) on 
"/.interval (0,1) 

epsilon = sqrt(((l-q~(2*L))/(l-q~(2*L+2)))*(q~2/(l+q~2))) ; 

°/ This is the value calculated 
°/ in Lemma 2.4.7 

gapest(i-9) = gamma * (l-sqrt(2)*epsilon) ~2; 

%This is the lower bound 
^provided in Theorem 2.4.4 

Deltavec(i-9) = Delta; 
H=[] ; 7 frees up memory- 
end 

plot (Deltavec , gap, 'b') 
hold on 

plot (Deltavec , gapest, 'r') 

We also include the programs Hmltn.m and SectrProj.m. 
Hmltn.m 

function Hmltn=Hmltn(L, J, Delta) 
%produces the spin J XXZ-Hamiltonian 



°/ A = 0.1; gives Delta = 1.005 



APPENDIX C. Matlab Code 



67 



Dellnv = l/Delta; °/„sqrt (1-A*A) ; 
A = sqrt(l-Delta~(-2)) ; 

y„Dimension of 1-site Hamiltonian 
diml = 2*J+1; 

%Pauli Spin Matrix in third (or Z) direction 
S3 = sparse (1: diml, 1: diml, J-(0:2*J) , diml, diml) ; 
%Spin-raising operator 
Splus = sparse (diml, diml) ; 
for j=0: (2*J-1) 

Splus = Splus + sparse(2*J-j ,2*J-j+l,sqrt((2*J-j)*(j+l)) , diml, diml) ; 
end 

%Spin lowering operator 
Sminus = transpose (Splus) ; 

°/oSS = Isotropic nearest neighbor interaction S(x).S(x+l) 

SS = kron(S3,S3) + (1/2) * (kron(Splus , Sminus) + kron(Sminus, Splus) ) ; 

°/ Ising = Ising Nearest neighbor interaction S3(x)*S3(x) 

Ising = kron(S3,S3); 

%h = nearest neighbor interaction 

h = -DelInv*SS - (1-Dellnv)*lsing+J~2*speye(diml~2) ; 

°/ H = sum over all nearest neighbor pairs of h(x,x+l) = Hamiltonian 

H = sparse(diml~L,diml~L) ; 

for x=l: (L-l) 

H = H + kron(eye(diml" (x-1) ) ,kron(h,eye(diml~(L-l-x)))) ; 
end 
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%Bdry = boundary term A*(S3(x+l)-S3(x)) 

Bdry = J*A* (kron(speye (diml" (L-l) ) ,S3) - kron(S3,speye (diml" (L-l)))) ; 
H = H +Bdry; 

Hmltn=H ; 

SectrProj.m 

function SectrProj=SectrProj (L, J, sector) 

°/„produces a projection matrix onto the the subspace of total S3 "sector" 
°/ for a spin J chain of length L 

downspins = J*L-sector; 

°/„Dimension of 1-site Hamiltonian 
diml = 2*J+1; 

°/„Pauli Spin Matrix in third (or Z) direction 
S3 = sparse (1: diml, 1: diml, J-(0:2*J) , diml, diml) ; 

%S3tot = total third-component of spin operator 
S3tot = sparse(diml"L,diml~L) ; 
for x=l:L, 

S3tot = S3tot + kron(speye (diml" (x-1) ) ,kron(S3 , speye (diml" (L-x) ))) ; 
end 
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°/ We now define the projection to the sector specified by 
% S3tot = ( J*L-downspins) 
Proj = speye(diml~L) ; 
for n=0: (2*J*L) 

if ne(n,downspins), 

Proj = Proj*(S3tot - (J*L-n)*speye(diml~L))/(n-downspins) ; 

end; 
end; 

°/„Proj is the orthogonal projection onto sector with specified number 

% of downspins . Next we want to define a projection from this subspace 

% to a vector space of the same dimension. 

°/„The command "find" finds the nonzero elements of Proj, 

°/ which is what we need. 

[11,12] = find(Proj) ; 

dim = length (II) ; 

NewProj = sparse (II , 1 : dim, ones (dim, 1) ,diml~L, dim) ; 



SectrProj=NewProj ; 
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